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ABSTRACT 

The Fermi gamma-ray satellite has recently detected gamma-ray emissions from radio galaxy cores. From 
these samples, we first examine the correlation between the luminosities at 5 GHz, L^cm, and at 0.1-10 GeV, 
L 7 , of these gamma-ray loud radio galaxies. We find that the correlation is significant with L 7 oc L 1 ^^ based on 
a partial correlation analysis. Using this correlation and the radio luminosity function (RLF) of radio galaxies, 
we further explore the contribution of gamma-ray loud radio galaxies to the unresolved extragalactic gamma- 
ray background (EGRB). The gamma-ray luminosity function is obtained by normalizing the RLF to reproduce 
the source count distribution of the Fermi gamma-ray loud radio galaxies. We find that gamma-ray loud radio 
galaxies will explain ~ 25% of the unresolved Fermi EGRB flux above 100 MeV and will also make a signif- 
icant contribution to the EGRB in the 1-30 MeV energy band. Since blazars explain 22% of the EGRB above 
100 MeV, radio loud active galactic nuclei (AGNs) populations explain ^47% of the unresolved EGRB. We 
further make an interpretation on the origin of the EGRB. The observed EGRB spectrum at 0.2-100 GeV does 
not show an absorption signature by the extragalactic background light. Thus, the dominant population of the 
origin of EGRB at very high energy (> 30GeV) might be nearby gamma-ray emitting sources or sources with 
very hard gamma-ray spectrum. 

Subject headings: cosmology: diffuse radiation - galaxies : active - gamma rays : theory 



1. INTRODUCTION 

The origin of the extragalactic diffuse MeV and GeV 
gamma-ray background (EGRB) radiation has been argued 
for a long time in astrophysics, although it is well known 
that radio quiet active galactic nuclei (AGNs) take account 
for the cosmic X-ray b ackground (CXB) below several hun- 
dred keV (for reviews lBoldtlll987t iFabian & Barconslll992l 
lUeda et alJIIOOl lHasinger et al.l l2005b iGilli et al.ll2007l) . The 
EGRB spectrum at 0.3-30 MeV is measured by SMM 
(IWatanabeet al.lll997h and C OMPTEL on board the Comp- 
ton Gamma-Ray Observatory (Kappada th et alj[l9 96). In the 
GeV energ y range, EGRB was first discovered by the SAS - 
2 satellite dFichtel et all Il978t iThompson & Fichtell 119821) . 
EGRET (Energetic Gamma-Ray Experiment Telescope) on 
board the Compton Gamma-Ray Ob servatory confirmed the 
EGRB spectrum at 0.03-50 GeV (ISreekumar et all 11998b 
Stron get al.l2004l) . Recently, LAT (Large Area Telescope) on 
board the Fermi Gamma-ray Space Telescope (Fermi) made 
a new measurement of the EGRB spectrum from 0.2 to 100 
GeV (lAbdo et alJ2010eh . The observed integrated EGRB flux 
(E > 100 MeV) is 1 .03 x 10" 5 photons/cm 2 /s/sr with a photon 
index of 2.41. 

Several sources have been suggested to explain the MeV 
background. One is the nuclear decay gamma-ray from 
Type la supernovae (SNe IalClavton & Warclll 1975b IZdziarskil 
1 19961: IWatanab e~et alJI 1999h . However, the recent measure- 
ment of the cosmic SNe la rate suggests that the expected 
flux from SNe la is one orde r of magn i tude lower than the 
measured MeV EGRB flux dAhn et alj 120051: IStrigari et ail 
12005b iHoriuchi & Beacomll2010h . Comptonization emission 
from non-thermal electrons in AGN coronae is also proposed 
dlnoue et al.ll2008l) . This model explains the origin of CXB 
and the MeV EGRB in the same population. Blazars, which 
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is one type of AGNs with the direction of a relati vistic jet co- 
incidi ng with our line o f sight, are also proposed (lAiello et ail 
120091) . Very recently, iMassaro & Ajellol (12011b has shown 
that the gamma-ray emission from lobes of radio galaxies 
will explain ~ 10% of the MeV background flux. MeV mass 
scale dark matter (DM) a nnihilations has also been discussed 
(lAhn & KomatsvJ l2005allbl) . but there is no natural particle 
physics candidate for a dark matter with a mass scale of MeV 
energies. However it is still discussing about the origin of the 
MeV background due to the difficulties of MeV gamma-ray 
measurements. 

In the case of the GeV background, sin ce blazars are 
domi n ant extragalactic g amma-ray sources (Hartm an et al.l 
11999b lAbdo et all 1201 Ogh . it is expected that unresolved 
population of blazars would explain the GeV EGRB 
dPadovani et alJl993blStecker et al.ll99llSalamon & Steckerl 
11994b IChiangetalJ 11995b IStecker & Salamonl [1996b 
Chi ang & Mukheried 11998b iMukheriee & Chiang! 
1999b iMiicke & Pohll 12000b iNarumoto & Totanil 12006b 
Giom mi et al.l 120061: iDermerl 120071: IPavlidou & V enters 
2008b lKneiske~& Mannhelrrl koOSl: IBhattacharva et alJl2009b 
llnoue & Totanill2009l) . Recently. lAbdo et alj d2010fl) showed 
that unresolved blazars can explain ~ 22% of the EGRB 
above 0.1 Ge V by analyzing 11 -month Fermi AGN catalog. 
Very recently IStecker & Ventersl (1201 Ol) proposed that the 
unresolved blazar population would be able to explain the 
EGRB spectrum below 1 GeV by taking into account of the 
energy dependence of source confusion effects. 

Other gamma-ray emitting extragalactic sources have also 
been discussed as the origin of the GeV EGRB. Those 
are intergalactic shocks produced by the large scale struc- 
ture formation (|Loeb & Waxmanl 120001: iTotani & Kitavamal 
l20MlMhiiatil2002l:lKeshet et alj2q03blGabici & Blasil2003l) 
normal and sta r burst g alaxies ( Pavlidou & Fi elds! 1 2002 
IThompson etail 12007b IBhattacharva & Sreekurriarf 12009 
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iMakiva et aLll20 1 Ob iFields et alj|2010h high Galactic latitude 
pulsa rs dFaucher-Giguere & Lo eb 2010: fSiegal-Gaskins et alJ 
120101) . kilo-parsec (kpc) size AGN je t s fstawarz et alT 2006), 
radio quiet AGNs (llnoue et all 120081; llnoue & Totaiiifl2009l) . 
and GeV mas s scale DM annih i lation or decay (s e e e.g . 
JungmanetaD 119961: iBergstroml ]2000t lUllio et al.1 12002b 



Oda et al. 2005 



Ando & Komatsu 2006; Horiuchi & Ando 



200 61: lAndoetalJ 120071; lAhnet alJ 120071: lAndol 120091: photon spectrum, dN/d 
Kawasaki et alj|2009l) . 

Here, Fermi has recently detected GeV gamma-ray emis- 
sions from 1 1 misaligned AGNs (i.e. radio galaxies), which 
are one type of AGNs with the directio n of a relativistic je t 
not coinciding with our line of sight (lAbdo et all [2010b). 
Their average photon index a t 0.1-10 GeV is ~ 2.4 which is 
same as that of Ge V EGRB dAbdo et all 1201 Pel) and blazars 
dAbdo et alJl2"0T03) . Although they are fainter than blazars, 
the expected number in the entire sky is much more than 
blazars. Then, it is naturally expected that they will make a 
significant contribution to EGRB. Therefore, in this paper, we 
study the contribution of such gamma-ray loud radio galaxies 
(not blazars) to EGRB. 

To study the EGRB contribution of gamma-ray loud radio 
galaxies, their gamma-ray luminosity function (GLF) is re- 
quired. Because of limited samples, it is not a straightfor- 
ward task to construct it using only Fermi samples. There- 
fore, we first investigate the correlation between radio and 
gamma-ray luminosities of gamma-ray loud radio galaxies. 
In the case of blazars, the correlation between gamma-ray 
and radio luminosities have been presented in many pa- 
pers s ince EGRET era (j Padovani et al .1 119931: IStecker et alJ 
19931: ISalamon & Steckerl Il994t iDondi & Ghisellinil 119951: 
Zhang et alJl200U iNarumoto & Totanil 120061: iGhirlanda et al.1 



11 -month catalog. Table. 1 lists the gamma-ray loud radio 
galaxy samples used in this study. It gives the object name, 
the First Source Catalog (1FGL) Femzz-LAT source name, 
redshift, gamma-ray photon flux above 0.1 GeV, photon in- 
dex at 0.1-10 GeV, radio flux at 5 GHz, spectral index at 5 
GHz, radio classification. The definition of photon index, T, 
and spectral index, a, is the index of power-law differential 



: , and T— 1, respectively. 
From Table. 1, the mean photon index of gamma-ray spec- 
tra at 0.1-10 GeV, T c , is 2.39 and t he spread is 0.28 This 
index is same as that of bla zars. 2.40 dAbdo et al.l l2010f). and 
the EGRB spectrum, 2.41 dAbdo et alJl2010eh. 

From the individual source studies dAbdo etal.1 l2009blcl 
1201 Ocb . the typical gamma-ray SEDs are well explained by 
synchrotron-self-Compton emission models. Here, by tak- 
ing account of particle cooling effect, the electron spectrum 
is given by dN /d^ e oc j~ p at j e < j\, T and dN/dj e °c 7~ (p+1) 
at 7 C > 7br, where 7b r is the cooling break lorentz factor. The 
inverse Compton (IC) photon spectrum is given by 



dN/de oc 



Hp+D/2 

-(p+2)/2 



£ < Cbr, 
e > Cbr, 



(1) 



2010aHr Since the radio luminos ity function (RLF) of ra- 
dio galaxies is wel l studied (see e.g. lDunlop & Peacock! 19901 : 
Will ott et alJ2001l) . we are able to obtain the GLF by convert- 
ing the RLF to the GLF using a luminosity correlation. With 
that GLF, we evaluate their contribution to EGRB. 

This work is organized as follows. Samples used in this 
study are shown in §f2] In §f3] we will investigate the ra- 
dio and gamma-ray luminosity correlation and determine our 
GLF to reproduce the number of detected gamma-ray loud 
radio galaxies. The EGRB will be calculated and compared 
with the observed data in §J4] Discussion and conclusion 
will be given in §[5] and §|6] respectively. Throughout this 
paper, we adopt the standard cosmological parameters of 
(ft,fij#,fi A ) = (0.7,0.3,0.7). 

2. SAMPLES 

Fermi has reported the detection of 1 1 Fanarof-Riley (FR) 
radio galaxies including 7 type-I FR (FRI) galaxies and 4 
type-II FR (FRII) galaxie s by the entire sky survey for 15 
months dAbdo et al.1 1201 Obi) . FRI galaxies have decelerating 
jets, knots at kpc distance from the core and edge-darkened 
lobes, while FRII galaxies have relativisti c jets and edge- 
brigh tened radio lobes with bright hotspots dFanaroff & Rileyl 
1974). In the scheme of the AGN jet unification scenario, FRI 
and FRII galaxies are the misaligned AGN populations of BL 
Lac o bjects and flat spectrum radio quasars (FSRQs), respec- 
tively dUrry & Padovanil 19951) . 

In this study, we use 10 samples (6 FRIs and 4 FRIIs) from 
lAbdo et ail d2010b|) which were al ready reported in 1 1 -month 
Fermi catalog dAbdo et al.ll2010 a g). This is because, as dis- 
cussed in §|3] b elow, we use the detection efficiency of Fermi 
shown in Abdo et al.l d2010j) which are constructed from the 



where eb r corresponds to the IC phot on energy from electrons 
with 7b r dRvbicki & Lightmanll 1979b . 

The SED fitting for NGC 1275 and M87 show that the IC 
peak ene rgy in the rest frame locates at ~ 5 MeV ( Ab do et ail 
2009b c). In this study, we use the mean photon index, T c , 
as r at 0.1-10 GeV and we set a peak energy, e^, in photon 
spectrum at 5 MeV for all gamma-ray loud radio galaxies as a 
baseline model. Then, we are able to define the average SED 
shape of gamma-ray loud radio galaxies for all luminosities 
as dN/de oc e" 2 39 at e >5 MeV, and dN/de oc e" 189 at e <5 
MeV by following Equation Q] 

However, only 3 sources are currently studied with multi- 
wavelength observational data. We need to make further stud- 
ies of individual gamma-ray loud radio galaxies to understand 
their SED properties in wide luminosity ranges. We examine 
other spectral models in the Section [ 



3. GAMMA-RAY LUMINOSITY FUNCTION 

3.1. Radio and Gamma-ray Luminosity Correlation 

To estimate the EGRB contribution of gamma-ray loud ra- 
dio galaxies, we need to construct a GLF. However, because 
of a small sample size, it is difficult to construct a GLF us- 
ing the current gamma-ray data only. Here, the RLF of ra- 
dio galaxies is extensively studied by previous works (see e.g. 
iDunlop & Peacock! 1 1990t I Willott etalJ 1200 ll) . If there is a 
correlation between the radio and gamma-ray luminosities, 
we are able to convert the RLF to the GLF with that cor- 
relation. In the case of blazars, it has been suggested that 
there is a correlation betw een radio and g amma -ray luminos- 
ity from the EGRET era dPadovani et al.ll 19931: fetecker et al.1 



1991 ISalamon & Steckerl 11994 IDondi & Ghisellinil 119951: 



Zhang et al. 2001; N arumoto & Tot ani 2006), although it has 



also been discussed that this correlation can not be firmly 



established because of flux limited samples (Mueckeetal 
119971). Recently, using the Fermi samples, IGhirlanda et al 



d2010allbl) confirmed that there is a correlation between the ra- 
dio and gamma-ray luminosities. 

To examine a luminosity correlation in gamma-ray loud ra- 
dio galaxies, we first derive the radio and gamma-ray lumi- 
nosity of gamma-ray loud radio galaxies as follows. Gamma- 
ray luminosities between the energies e\ and e^ are calculated 
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TABLE 1 

Observed parameters of gamma-ray loud radio galaxies. 



Object Name 


1FGL Name 


z 

[xlO- 9 


F-y 

ph/cm 2 /s] 


r 


5r 

[Jy] 




Class 


REF 


3C78/NGC 1218 


1FGLJ 0308.3+0403 


0.029 


4.7±1.8 


1.95±0.14 


0.964 ± 0.048 


0.64 


FRI 


1 


3C 84/NGC 1275 


1FGLJ 0319.7+4130 


0.018 


222 ±8 


2.13 ±0.02 


3.10±0.02 


0.78 


FRI 


2 


3C 111 


1FGLJ 0419.0+3811 


0.049 


40±8 


2.54±0.19 


1.14 ±0.0° 


-0.146 


FRII 


3 


PKS 0625-354 


1FGLJ 0627.3-3530 


0.055 


4.8± 1.1 


2.06±0.16 


0.60 ±0.03 


0.53 


FRI 


1 


3C207 


1FGLJ 0840.8+1310 


0.681 


24±4 


2.42 ±0.10 


0.51 ±0.02 


0.9 


FRII 


2 


PKS 0943-76 


1FGLJ 0940.2-7605 


0.27 


55 ±12 


2.83±0.16 


0.79 ±0.03 


0.79 


FRII 


4 


M87/3C 274 


1FGLJ 1230.8+1223 


0.004 


24±6 


2.21 ±0.14 


4.0 ±0.04 


0.79 


FRI 


2 


Cen A 


1FGLJ 1325.6-4300 


0.0009 


214±12 


2.75 ±0.04 


6.98 ±0.21 


1.2 


FRI 


1 


NGC6251 


1FGLJ 1635.4+8228 


0.024 


36±8 


2.52±0.12 


0.35 ±0.045 


0.72 


FRI 


2 


3C 380 


1FGLJ 1829.8+4845 


0.0692 


31 ± 18 


2.51 ±0.30 


7.45 ±0.047 


0.71 


FRII 


2 


NOTE. — 1FGL Name: the First Source Catalog (IFGL) Fermi-LAT 


source name 


z: redshift of the 


source, F~ f : gamma-ray photon flux above 100 



MeV in 10 photons/cm /s, Y: photon index at 0.1-10 GeV, Sr: radio flux densit y at 5 GHz in Jy , a r : radio spectral i ndex at 5 GHz, Clas s : FRI is type 
I of Fanaroff-Riley galaxy and FRII is type I I of Fanaroff-Riley galaxy i Fanaroff & Riley 1974). Refer ences — ( l)IUnger et al 1 1984); Saikia^^l 
i 1986); Baum et al. i 1988); Ekers et al. 1 1989); Jones & McAdam i 1992); Burns et al. 1 1983); Moreanti et al. < 1993); (2) Bennetl 1962); Spinrad et al 



1985); Laing et al. 1 1983); (3) Linfield & Perley I 
1 No error is reported in Linfield & Perley ( 1984). 



by 



L 7 (ei 



e 2 ) = 4ird L {z) 



2 S 7 (ei,e 2 ) 



(1+z) 



2-r ' 



; (4) Burgess & Hunstead i 2006a b). 



(2) 



where c/l(z) is the luminosity distance at redshift, z, T is the 
photon index and ^(ei , e 2 ) is the observed energy flux between 
the energies ei and e 2 . The energy flux is given from the pho- 
ton flux Ey, which is in the unit of photons/cm 2 /s, above t\ 
by 



(r-l)ei 

2-r 



2-r 



■1 



5 7 (ei,e 2 ) = eiln(e 2 /ei)F 7 , (r = 2) 



f 7 , cry 2) (3) 



(4) 



Radio luminosity is also calculated in the same manner. 

Figure Q] shows the 5 GHz and 0.1-10 GeV luminosity re- 
lation of Fermi gamma-ray loud radio galaxies. Square and 
triangle data represents FRI and FRII radio galaxies, respec- 
tively. The solid line shows the fitting line to all the data. The 
function is given by 

log 10 (L 7 ) = (-3.90 ± 0.61) + (1 . 16 ± 0.02) log 10 (L 5GHz ), (5) 

where errors show l-a uncertainties. In the case of blazars, 
the slope of the correlation between L 7 (> lOOMeV), lumi- 
nosity above 100 MeV, and radio luminosity at 20 GHz is 
1.07 ±0.05 dGhirlanda et al.ll2010bl) . The correlation slopes 
of gamma-ray loud radio galaxies are similar to that of 
blazars. This may indicate that emission mechanism is similar 
in gamma-ray loud radio galaxies and blazars. 

We need to examine whether the correlation between the 
radio and gamma-ray luminosities is true or not. In the flux 
limited observations, the luminosities of samples are strongly 
correlated with redshifts. This might result in a spurious lu- 
minosity correla t ion. As in previous works on blazar sample s 
(lPadovanilll99l IZhang etalJl200lt iGhirlanda et all 1201 Obi) , 
we perform a partial correlation analysis to test the correla- 
tion between the radio and gamma-ray luminosities exclud- 
ing the redshift dependence (see the Appendix for details). 
First, we calculat e the Spearman r ank-order correlation co - 
efficients (see e.g. iPress et al.llT992b . The correlation coeffi- 
cient is 0.993, 0.993, 0.979 between log 10 L 5G Hz and log 10 L 7 , 

and 



between log 10 L 5 GHz an d redshift, and between log 10 L 7 
redshift, respectively. Then, the partial correlation coefficient 
becomes 0.866 with chance probability 1.65 x 10~ 6 . There- 
fore, we conclude that there is a correlation between the radio 
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FIG. 1 . — Gamma-ray luminosity at 0.1-10 GeV versus radio luminosity 
at 5 GHz. The square and triangle data represents FRI and FRII galaxies, 
respectively. The solid line is the fit to all sources. 

and gamma-ray luminosities of gamma-ray loud radio galax- 
ies. 

3.2. Gamma-ray Luminosity Function 

In this section, we derive the GLF of gamma-ray loud radio 
galaxies, p 7 (L 7 ,z)- There is a correlation between the radio 
and gamma-ray luminosities as Equation|5] With this correla- 
tion, we develop the GLF by using the RLF of radio galaxies, 
p r (L r ,z), with radio luminosity, L r . The GLF is given as 



■ K——p r (L r ,z), 
aL-y 



(6) 



where k is a norma lization factor. We use the 151 MHz RLF 
dWillott et al.ll200l . Since they presented the formula of FRI 
and FRII RLFs separately, we combined them as in their pa- 
per because it is difficult to analyze each population separately 
with our limited number of samples. Moreov er, since the cos- 
mological parameters in Will ott et al.l (12001 1) are Qm = = 
and h = 0.5, we also convert the RLF to the standard cosmol- 
ogy adopted in this study. The RLF is given by 

p r (L r ,z) = r](z) x [p r ,mi(L r ,z)+ pr,mn(L r ,z)], (7) 

where r/(z) is the conversion factor of cosmology, /O r FRi is the 
FRI RLF, and p r FRII is the FRII RLF. The FRI RLF and FRII 
RLF are given by 



Pr,FRl(L r ,z) = 



.^°(fc) ai ^p(§)(l+Zl,c) fe Z>Zl,, 



(8) 
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TABLE 2 
The parameters of the RLF 



Willottetal. (2001) 



lo glo(Pu)) 

l°gioCZ&) 

Zl,c 
*1 

lo gio(Pn,o) 
«n 

%io(Lf Ic ) 
zir,c 
zn,i 

ZII.2 



-7.523 
0.586 
26.48 
0.710 

3.48 
-6.757 

2.42 
27.39 

2.03 
0.568 
0.956 



1 In units of Mpc 3 . 
' In units of W/Hz/Sr. 



exp 



^II.c 



Here, fu(z) is the evolution function and given by 



/n(z). (9) 



Mz)-- 



(10) 



T he parameters of these RLFs are summarized in Table. |2] As 
in lStawarz et al.l (120061) . the conversion factor of the cosmol- 
ogy r?(z) is 

_ d 2 V w /dndz 



Viz) 



(11) 



d 2 V/dVLdz ' 

The co moving volume element of cosmology in lWillott et alj 
J2001 is 

c i z 2 (2+z) 2 



d 2 V w 



dVLdz ml w (\+zf 



(12) 



where c is the speed of light and Hqw is 50km/s/Mpc. That 
of our standard cosmology is 



d 2 V 



cddzf 



dndz H Q (i+z) 2 V(i-n M -n A )(i+z) 2 +n M (i+z) 3 +n A : 

(13) 

where Ho is 70km/s/Mpc. 

Since Equation [5] is for radio luminosities at 5GHz in the 
unit of erg/s, we assume spectral index a r = 0.8 for all radio 
galaxies to co nvert 151 MHz lum inosity to 5GHz luminosity 
as assumed in lWillott et al.l (120011) . Although a r would affect 
the fraction of gamma-ray loud radio galaxies in radio galaxy 
population, a r does not affect main results on the EGRB cal- 
culation in this paper because our GLF is normalized to the 
cumulative source count distribution of the gamma-ray loud 
radio galaxies detected by Fermi. 

3.3. Source Count Distribution 

The normalization factor At, which corresponds to the frac- 
tion of gamma-ray loud radio galaxies against all radio galax- 
ies, is determined by the normalizing our GLF to the source 
count distribution of the Fermi radio galaxies, which is some- 
times called logN-logS plot or cumulative flux distribution. 
Source count distribution is calculated by 



N(> F 7 ) = 4tt 



dz 



d 2 V 
dildz 



t/L 7 p 7 (L 7 ,z), (14) 
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FIG. 2. — The source count distribution of gamma-ray loud radio galaxies 
in the entire sky. Solid and dashed curve corresponds to the all radio galaxies 
and gamma-ray loud radio galaxies, respectively. The data points shows the 
Fermi data after the conversion of detection efficiency. The error bar shows 
l—cr statistical uncertainty. 



where L 7 (z,F 7 ) is the gamma-ray luminosity of a blazar at 
redshift z whose photon flux at > 100 MeV is F 7 . Hereinafter, 
we assume z max = 5 and L 7 , max = 10 48 erg/s in this study. These 
assumptions hardly affect the results in this study. 

Since the completeness of the Fermi sky survey depends on 
the photon flux and photon index of a source, we need to take 
into account this effect (so called the detection efficiency) to 
compare GLF with the cumulative source count distribution 
of gamma-ray loud radio galaxies. The dete c tion ef ficiency 
of Fermi is shown in Figure 7 of lAbdo et all (l2010fj) for the 
sources in the 11-month catalog with test statistics TS> 50, 
at the Galactic latitude \b\ > 20°, and with a mean photon in- 
dex of 2.40, which is similar to that of gamma-ray loud radio 
galaxies (see Section |2]). It is shown that results for blazar 
source count distribution analysis did not change even if sam- 
ples with \b\ > 15° are included. Furthermore, even if they in- 
clude samples with TS> 25, the systematic uncertainties are 
small. Therefore, w e adopt the detection efficiency shown in 
lAbdo et ail d2010ii) to our samples from the Fermi 1 1 -month 
catalog, although not all the samples locates \b\ > 20° or have 
TS> 50. 

Figure|2]shows the source count distribution of gamma-ray 
loud radio galaxies. The data is after the conversion of the 
detection efficiency. Solid line shows in the case of At = 1 
which corresponds to the case that all radio galaxies emit 
gamma-rays. Dashed curve corresponds to the GLF fitted 
to the Fermi data with re = 0.081 ±0.011. — 1000 gamma- 
ray loud radio galaxies are expected with 100% complete en- 
tire sky survey above the flux threshold F 7 (> lOOMeV) = 
1.0 x 10" 9 photons cm" 2 s" 1 above 100 MeV. We note that 
the current detection efficiency of Fermi at F 7 (> lOOMeV) = 
1.0 x 10" 9 photons cm" 2 s" 1 is - 10" 3 . 

4. EXTRAGALACTIC GAMMA-RAY BACKGROUND 

We calculate the EGRB spectrum by integrating our GLF 
in the redshift and luminosity space, using the SED model 
shown in Section. [2] The EGRB spectrum is calculated as 





dt 




/ dz 




J 


lo 


dz 


Jl 



d 2 F(e)_ c 

dedVL ~ An 

dL[L 1 ,(\+z)e] 

X dt 
xexp[-r 7i7 (e,z)], 



x{1.0-o;(F 7 [L 7 ,z])} 



(15) 
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where t is the cosmic time and dt/dz can be calculated by 
the Friedmann equation in the standard cosmology. The min- 
imum gamma-ray luminosity is set at L i m ; m = 10 erg/s be- 
cause there is no reported gamma-ray loud radio galaxies be- 
low this value. w(F 7 [L 7 ,z]) is the detection efficiency of 
Fermi at the photon flux F 1 which corresponds to the flux 
from a source with a gamma-ray luminosity L 7 at redshift z. 

High energy 7-rays (>20 GeV) propagating the universe 
are absorbed by the interaction with the extragalactic back- 
ground light (E BL), also called as cosmic optical and infrared 
background. (ISalamon & Steckerl 119981; iTotani & Takeuchi 



2002; Kneiske et alJl2004; Stecker et al. 2006; Mazin & Rauel 



20071: iRaue & Mazinl 120081: iF ranceschi ni et al. 



Gilmore et a l. 2009; Fink e et aU 



2008 



201 



Razzaque et al. 200' 
Kneiske & Dole 2010). r 7i7 (e,z) is the optical depth of this 
background radiat ion. In this study, we adopt the model of 
Fink e etal.1 (1201 Ol) for EBL and r 7 a . 

The gamma-ray absorption creates electron-positron pairs . 
These pairs scatter the cosmic microwave background ra- 
diation to make the secondary emission component (so 
called cascade emission) to the absorbed primary emis- 
sion (Aharonia n et alj 119941: IWang et al.l 120011; iDaietalJ 
2002t iRazzaaue et al.l | 2004l: lAndoll2004t iMurase et al.ll2007b 
Kneiske & Man nheiml 120081 : llnoue & Totanj 120091; IVenteri 
20101) . We take into account the first generation of 



the cascade emission follo wing the formulations as in 
Kneiske & Mannheim ( 2008). In the following result, the cas- 
cade emission cause a small effects on the EGRB flux. Hence, 
the other generations of cascade emission do not have a seri- 
ous effects on our conclusion in this study. 

Fig. [3] shows the vl v EGRB spectrum in the unit of 
MeV 2 cm" 2 s" 1 MeV" 1 sr" 1 predicted by our GLFs. The in- 
trinsic (the spectrum without the EBL absorption effect), 
absorbed, and cascade components of the EGRB spec- 
trum and the total EGRB sp ectrum (absorbed+c ascade) are 
shown. The data o f HEA O-I dGniber et al.lll999h . Swift-BAT 



(Aiello et al. 2008), SMM (Watanabe et al. 1997), COMPTEL 



Kappadathet alJI 19961) . and Fenni-LAT (lAbdo et al.ll2010d) 



are also shown. As in the figure, the cascade emission does 
not contribute to the EGRB spectrum significantly. 

The expected EGRB photon flux above 100 MeV 
from gamma-ray loud radio galaxy populations is 0.26 x 
10~ 5 photons cm" 2 s" 1 sr" 1 . As the unresolved Fermi EGRB 
flux above 100 M eV is 1.03 x 10~ 5 photons cm" 2 s" 1 sr" 1 
dAbdo et all 1201 Pet) , the gamma-ray loud radio galaxies ex- 
plains ^25% of the unresolved EGRB flux. For the com- 
parison, recent analysis of Fermi blazars showed that blazar s 
explains ~ 22% of the unresolved EGRB (lAbdo et al.ll2010fl) . 
Therefore, we are able to explain ~ 47% of EGRB by radio 
loud AGN populations. 

To avoid the instrument dependence, we also evaluate the 
total EGRB photon flux (i.e. resolved + unresolved) from 
gamma-ray loud radio galaxies. The contribution to total 
EGRB is 0.27 x 10" 5 photons cm" 2 s" 1 sr" 1 and this corre- 
sponds to ~ 19% of total Fermi EGRB flux w hich is 1.42 x 
10" 5 photons cm" 2 s" 1 sr" 1 dAbdo et al.ll2010el) . 

Figure [3] also shows that gamma-ray loud radio galaxies 
would also contribute to the MeV EGRB at > 1 MeV signifi- 
cantly. Although it has been sugges ted that radio quie t AGNs 
( llnoue et all 120081) . MeV blazars (lAiello et all 120091). MeV 
DM annihilation (e.g. lAhn & Komatsull2005al) . it is still un- 
certain because of luck of observational evidences. As shown 
in this study, gamma-ray loud radio galaxies would also be a 
candidate for the origin of MeV background. This situation 



will be solved by future X -ray and MeV gamma -ray experi- 
ments such as ASTR O-H 1 dTakahashi et al.ll2010h and DUAL 
gamma-ray mission (Bo ggs et al.l l2010). respectively. 

We should examine the uncertainties in the model predic- 
tion. Since the normalization of the GLF is determined from 
10 samples, there is a statistical uncertainty of 32% in its nor- 
malization of the EGRB at 68 % confidence level. The corre- 
lation of the radio and gamma-ray luminosities also has uncer- 
tainties in their slope and normalization as in § 13.11 By taking 
into those uncertainties, the fraction of gamma-ray loud radio 
galaxies in unresolved EGRB varies from ~ 10% to ~ 63%. 
Furthermore, as discussed in Stecker & Venters (2010|), the 
energy dependence of source confusion effects would alter the 
EGRB spectrum below ^1 GeV. However, the angular reso- 
lution is also depends on the p osition of the source in the field 
of view (lAtwood et al.ll2009l) . The further careful evaluation 
is required to discuss the source confusion effects. 

5. DISCUSSION 
5.1. Fraction of Gamma-ray Loud Radio Galaxies 

Since jet is brighter to observers with a smaller viewing an- 
gle from the jet axis because of the beaming effect, the frac- 
tion of gamma-ray loud radio galaxies k would be related with 
the vie wing angle. It is believe d that radio galaxies have bipo- 
lar jets (lUrry & Padovanill 1995b . The fraction of radio galax- 
ies with the viewing angle < 9 is given as k = (1 -cosf?). In 
this study, the fraction of gamma-ray loud radio galaxies is 
derived as k = 0.081 as discussed in § 13.31 Then, the expected 
9 is < 24°. The viewing angle of NGC 1275, M 87, and Cen 
A is derived as 2 5°, 10°, and 30° by SED fitting dAbdo etalj 
I2009bld 120 lOd) . respectively. Therefore, our estimation is 
consistent with the observed results. 

Here, beaming factor <5 is defined as r _1 (l - /3cos60 _1 , 
where T is the bulk lorentz factor of the jet and (3 = yl — 1 /T 1 . 
If T ~ 10 which is typical for blazars, 5 becomes ~1 with 
9 = 24°. This value means no significant beaming effect be- 
cause the observed luminosity is <5 4 times brighter than that in 
the jet rest frame. On the other hand, if 2 < T < 4, <5 becomes 
greater t han 2 with 9 = 24° (i.e . beaming effect becomes im- 
portant). Ghisellini et al. (2005) proposed the spine and layer 
jet emission model, in which jet is composed by a slow jet 
layer and a fast jet spine. The difference of T between blazars 
and gamma-ray loud radio galaxies would be interpreted with 
structured jet emission model. 

We note that k depends on a r as in Sec. 13.21 By changing 
a r by 0.1 (i.e. to 0.7 or 0.9), n and 9 changes a factor of 1.4 
and 1.2, respectively. Thus, even if we change a r , beaming 
effect is not effective with T ~ 10 but with lower T value, 

2<r<4. 

5.2. Uncertainty in the Spectral Modeling 

As pointed in Section [2] there are uncertainties in SED 
modeling because of small samples such as photon index 
(r) and break photon energy (eb r ). In the case of b lazars, 
IStecker & Salamonl (1 19961) ; (pavlidou & Ventersl (120081) calcu- 
lated the blazar EGRB spectrum including the distribution 
of the photon index by assuming Gaussian distributions even 
with ~ 50 samples. We performed the Kolomogorov-Smirnov 
(K-S) test to see the goodness of t he gaussian fit to our sampl e 
and check whether the met hod in lStecker & Salamonl (119961) : 
iPavlidou & Ventersl d2008l) is applicable to our sample. The 

1 ASTRO-H : http://astro-h.isas.jaxa.jp/index.html.en 
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Photon Energy [MeV] 

FIG. 3 . — EGRB spectrum from gamma-ray loud radio galaxies in the unit of MeV 2 cirT 2 s~' MeV~' sr~' . Dashed, dotted, dot-dashed, and solid cu rves show the 
intrin sic spectrum (no absorption), ab sorbed, cascade, and total (ab sorbed+cascade) EGRB spectrum, re spectively. The o bserved data of HE AO-I ( Gr uber et al.l 
GUI), Swift-BAT jAiello et alj|2 008), SMM jWatanabe et"aT|[l997h . COMPTEL I Kappadath et al. 1996), and Fenni-LAT JAbdo et aU20l0el) are also shown by 
the symbols indicated in the figure. 



chance probability is 12%. This means that the Gaussian dis- 
tribution does not agree with the data. To investigate the dis- 
tribution of photon index, more sample would be required. 

We evaluate the uncertainties in SED models by using var- 
ious SEDs. Figure [4] shows the total EGRB spectrum (ab- 
sorbed + cascade) from the gamma-ray loud radio galaxies 
with various photon index and break energy parameters. The 
contribution to the unresolved Fermi EGRB photon flux above 
100 MeV becomes 25.4%, 25.4%, 23.8% for T = 2.39, 2.11, 
and 2.67, respectively. In the case of V = 2.11, the contri- 
bution to the EGRB flux above 10 GeV becomes significant. 
For the MeV background below 10 MeV, the position of the 
break energy and the photon index is crucial to determine the 
contribution of the gamma-ray loud radio galaxies. As shown 
in the Figure |4] higher break energy and softer photon index 
result less contribution to the MeV background radiation. To 
make a further discussion on the SED modeling, the multi- 
wavelength spectral analysis of all GeV observed gamma-ray 
loud radio galaxies are required. 

5.3. Flaring Activity 

It is well k nown that blazars are var iable sources in gamma- 
ray (see e.g. HjxEetaDEooH Hojol). If gamma-ray loud 
radio galaxies are the misalign ed population of blaz ars, they 
will also be variable sources. Kataoka et al.1 (1201 Oh has re- 
cently reported that NGC 1275 showed a factor of ~ 2 vari- 
ation of gamma-ray flux. For other gamma-ray loud radio 
galaxies, such a si gnificant variation has not observed yet 
( Abdo ^taTI 1201 Ob i). Therefore, it is not straightforward to 
model the variability of radio galaxies currently. In this pa- 




10" 2 10" 1 10° 10 1 10 2 10 2 10 4 10 5 10 6 

Photon Energy [MeV] 



FIG. 4. — Gamma-ray loud radio galaxy EGRB spectra in the unit of 
[MeV 2 cm~ 2 s MeV~' sr~'] for various SED parameters, T and £(,r- The 
curves shown in this figure is the total EGRB spectrum (absorbed + cascade). 
Solid, dashed, dotted, dot-dashed, double-dotted dashes curves corresponds 
to the EGRB spectrum with (T, e br ) =(2.39, 5.0 MeV), (2. 1 1 , 5.0 MeV), (2.67, 
5.0 MeV), (2.39, 0.5 MeV), and (2.39, 50.0 MeV), respectively. The observed 
data shown here are same as Figure [3] 



per, we used the time-averaged gamma-ray flux of gamma-ray 
loud radio galaxies in the Fermi catalog, which is the mean 
of the Fermi 1-year observation. More observational informa- 
tion (e.g. frequency) is required to model the gamma-ray vari- 
ability of radio galaxies. Further long term ferm/observation 
will be useful and future ground based imaging atmospheric 
Cherenkov Telescope, Cherenkov Telescope Array (CTA) 2 

2 CTA: http://www.cta-observatory.org/ 
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would be a key to understanding short period variabilities. 

5.4. Origin of the GeVEGRB 

In this study, we find that the contribution of gamma-ray 
loud radio g alaxies to the unreso lved EGRB above 100 MeV 
is ^25 %. lAbdo et all d2010fl) recently showed that unre- 
solved blazars can explain only ~ 22 % of the unresolved 
EGRB by analyzing one year catalog of the Fermi blazars. 
Therefore, the origin of rest ^53% of EGRB is still missing. 

Various gamma-ray emitting extragalactic sources have 
also been discussed as the GeV EGRB origin. Those are 
intergalactic shocks produced by the l arge scale structure 
forma t ion dLoeb & Waxmanl 120001: iTotani & Kitavamal 
2000t iMiniatil l2002t iKeshetetalJ (20031: iGabici & Blasl 



Total 
z = 0.0-0.8 
z = 0.8-2.0 
z = 2.0-5.0 
Fermi-LAT 



2003 
2007 
2010 



galaxies ( Pavlidou & Fields 2002; Thompson et al. 
iBhattach arva & Sreekumar 120091 iMakiva et al.l 
iFields et all 120101; IStecker & Ventersl BO 10b high 
Galactic latitude pulsars (IFaucher-Giguer e & Loebl 1201 Ol 
ISiegal-Gaskins et al.ll20Tob , kilo -parsec (k pc) size AGN jets 
dStawarz et al.l 120061) . radio quiet AGNs (llnoue et all 120081; 



or decay (see e.g. Jungman et al. 


1996; Bergstrom 
Ando & Komatsu 


12000 


Ullio et all 2002; 


Odaetal. 2005; 


2006 


Horiuchi & Ando 


2006; Ando et al. 2007; Ahn et al. 


2007 


Ando 2009; Kawasaki et al. 2009). 



Fig. [5] shows the gamma-ray loud radio galaxy EGRB 
spectra at each redshift bins. Because of EBL, the spectrum 
above > 30 GeV shows the absorbed signature. Here, the 
cosmological sources have their evolution peaks at z = 1 ~ 2 
suc h as cosmic star formation history and AGN activity (see 
e.g. lHopkins & B eacom 2006; Ued a et al]|2003b . This means 
that the gamma-rays from extragalactic sources (e.g. galax- 
ies and AGNs) will experience the EBL absorption. How- 
ever, as shown in Fig. [5] Fermi EGRB spectrum does not 
show such an absorbed signature. This might suggest that 
nearby gamma-ray emitting sources or sources with very hard 
gamma-ray spectrum would be the dominant population of 
EGRB above 10 GeV. To address this issue, we should await 
the EGRB information above 100 GeV by future observations 
such as Fermi. CTA would also be possible to see the EGRB 
at much higher energy band. We also need to examine the 
EBL models at high redshift. It is expected that CTA will 
see blazars up to z ~ 1 .2 at very high energy band > 30 GeV 
( llnoue et al.ll2010h . Therefore, Fermi and CTA will be a key 
to understanding the origin of EGRB. 

5.5. Implication to the AGN Unification Scenario 

AGN unification scenario exp lains various propertie s of 
AGNs in terms of viewing ang le (lUrrv & Pado vani 1995b . In 
the scheme of AGN jet unification scenario, FRI and FRII 
galaxies are thought as misaligned populations of BL Lacs 
and FSRQs, respectively. 

From Table. Q] the mean photon index of FR Is and FRIIs 
are 2.27 and 2.58, respectively. Therefore, the FRI popula- 
tion tend s to have ha r der spe ctra than the FRII population as 
shown in I Abdo et al.l d2010bb . This trend is also same as that 
between BL Lacs and FSRQs dAbdo et alj|2010gh . This result 
would support that FRIs and FRIIs are the misaligned popu- 
lation of BL Lacs and FSRQs. 

It is also important to compare the cosmological evolutions 
of blazars and radio galaxies based on the re cent Fermi data. 
Altho ugh a theoretical blazar GLF model (In oue & Totanil 
2009) is briefly compared with the Fermi EGRB result and the 
cumulative source count distribution of the Fermi blazar data 




10* io- 
Photon Energy [MeV] 

FIG. 5. — Gamma-ray loud radio galaxy EGRB spectra in the unit of 
[MeV 2 cm~ 2 s MeV~' sr _1 ] at each redshift bins. Solid, dashed, dotted, 
and dot-dashed curves shows the EGRB spectrum at z = 0.0-5.0, z = 0.0—0.8, 
z — 0.8-2.0, and z — 2.0-5.0, respectively. The thin and thick curve corre- 
sponds to the intrinsic spectrum (not absorbed) and the total sp ectrum (ab- 
sorbed + cascade), respectively. The observed data of Fermi-LAT I A bdo et all 
2010e) is also shown. 



dlnoue etal.|[2010t llnoue & Totanil 1201 U llnoue etalJl20TTh . 
comparison in redshift space has not yet been done. This is 
because redshif ts of about a half o f BL Lac samples has not 
yet determined dAbdo et al.ll2010g|). 

In addition, llnou e & Totani ( 2009) treated blazars as 
one population u sing blazar sequence dFossati et al.l [1998; 
iKubo et al]|"l998l) . Blazar GLF models which di vide FSRQs 
and BL Lacs is required such as Dermei] (12007b to interpret 
the unification scenario. Since our model in this paper also 
does not treat FRI and FRII separately because of small sam- 
ples, GLF models of gamma-ray loud radio galaxies dividing 
these two population are also required. 

Therefore, redshift information of all blazars and more data 
of gamma-ray loud radio galaxies would be required to make 
a comparison of the cosmological evolutions of blazars and 
radio galaxies. 

6. CONCLUSION 

In this paper, we studied the contribution of gamma-ray 
loud radio galaxies to the EGRB by constructing their GLF. 
First, we explored the correlation between the radio and 
gamma-ray luminosities of gamma-ra y loud radio galaxie s 
which are recently reported by Fermi dAbdo et aT]|2"0 10g b). 
We found that there is a correlation L 7 oc L^'q^ 002 by a partial 
correlation analysis, where L 7 is the 0.1-10 GeV gamma-ray 
luminosity and Lsghz is the radio luminosity at 5 GHz. This 
slope index is similar to that of blazars. 

Based on this correlation, we defined the GLF of gamma- 
ray loud radio galaxies using the RLF of radio galaxies. We 
normalized the GLF to fit to the cumulative flux distribu- 
tion of Fermi sam ples by using Fermi detection efficiency 
dAbdo et alj|2010li) . Then, we predicted the contribution of 
gamma-ray loud radio galaxies to the MeV and GeV EGRB. 
The absorption by EBL and the reprocessed cascade emission 
are also taken into account. We found that gamma-ray loud 
radio galaxies will explain ^25% of the EGRB flux above 
100 MeV and also make a significant contribution to the 1-30 
MeV EGRB. Since blazars explain ~ 22% of EGRB, we are 
able to explain ~ 47% of EGRB by blazars and gamma-ray 
loud radio galaxies. 

We also make an interpretation on the origin of the EGRB 
above 10 GeV from the point of view of the EBL absorp- 
tion effect. Since the EBL absorption signature is still not 
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appeared in the EGRB spectrum, the origin would be nearby 
sources or sources with hard gamma-ray spectrum. We should 
await the EGRB data at higher energy band for this issue. 
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APPENDIX 
PARTIAL CORRELATION ANALYSIS 

To study the correlation between luminosities at different wavelengths, we use luminosities directly. However, the correlation 
in luminosity space is distorted by redshift, if samples are flux-limited. This will result in a spurious correlation. Therefore, we 
need to test t he correlation ex cluding the redshift dependence. Partial correlation analysis is the analyzing method for such a 
condition (see Padovani 1992, for details). The partial correlation analysis method is as follows. 

We have parameter sets of {x, • , y, z,), i = 1 , 2 . . . . , N. Let X, be the rank of jc, among the other x's, Yi be the rank of y, among the 
other y's, and Z, be the rank of z, among the other z's. The Spearman rank-order correlation coefficient between x and y is defined 
to be the linear correlation coefficient of the ranks as 

, , (A1) 



where X and Y is the mean of the X and the Y , respectively (see lPress et al.ll992l for details). The correlation coefficients between 
x and z, and y and z are also given in a same way. Then, the correlation coefficient between x and y excluding the dependence on 
the third parameter of z is evaluated as 

r w = rxy ~ Vyz , (A2) 
where r xz and r yz is the correlation coefficient between x and z and between y and z, respectively (Ke ndall & Stuarl| [T979). 



